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ABSTRACT 

Any viable theory of modified Newtonian dynamics (MOND) as modified gravity is 
likely to require fields in addition to the usual tensor field of General Relativity. For 
these theories the MOND phenomenology emerges as an effective fifth force probably 
associated with a scalar field. Here I consider the constraints imposed upon such theo- 
ries by solar system phenomenology, primarily by the absence of significant deviations 
from inverse square attraction in the inner solar system as well as detectable local pre- 
ferred frame effects. The current examples of multi-field theories can be constructed 
to satisfy these constraints and such theories lead inevitably to an anomalous non 
inverse-square force in the outer solar system. 
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1 INTRODUCTION 



In modified Newtonian dynamics (MOND) it is postulated 
that the effective gravitational acceleration, g, deviates from 
the Newtonian value, qm, below a critical acceleration, ao, 
in the sense that g « y/a^gN (Milgrom 1983). That the em- 
pirically determined value of ao (~ 10~* cm/s^) coincides 
with cHq to within an order of magnitude was immediately 
noticed by Milgrom who speculated that MOND may re- 
flect the influence of cosmology on local particle dynamics. 
In the context of General Relativity (GR), there is no such 
cosmological influence. This is essentially due to fact that 
GR embodies the equivalence principle in its strong form 
which forbids environmental influence on local dynamics, 
apart from tides. However, if the gravitational force is par- 
tially mediated by a long range scalar fleld, as, for example, 
in Brans-Dicke theory, it is no longer the case a local sys- 
tem is immune from cosmological influence. The scalar field, 
determined by the universal mass distribution and its time 
evolution, pervades the Universe and influences the gravita- 
tional dynamics of every subsystem. In Brans-Dicke theory 
this influence is evidenced by the cosmic evolution of the 
effective gravitational constant (Brans & Dicke 1961). 

This suggests that MOND may have its basis in scalar- 
tensor theory; indeed, the first relativistic theories proposed 
for MOND were scalar-tensor theories with non-standard as- 
pects: the aquadratic Lagrangian theory, AQUAL (Beken- 
stein & Milgrom 1984) , in which the scalar field Lagrangian 
is a general function of the usual scalar field invariant 
{F{(f}'°'(f}^a)), and phase-coupling gravity, or PCG, (Beken- 
stein 1988) in which the scalar fleld is complex with stan- 
dard Lagrangian but only phase coupling to matter. Both 



of these early attempts contain pathologies- superluminal 
propagation or instability of the background (Bekenstein & 
Milgrom 1984, Bekenstein 1990) . Moreover, in these theories 
the scalar field is assumed to couple to matter jointly with 
the gravitational, or Einstein, metric in order to preserve 
the universality of free fall (Weak Equivalence Principle or 
WEP). But if that coupling is conformal, as in Brans-Dicke 
theory, then there is no enhanced gravitational defiection of 
photons due to the scalar field. This is in dramatic confiict 
with observations of lensing by clusters of galaxies (Beken- 
stein & Sanders 1994). 

The lensing contradiction led to the idea that the re- 
lation between the Einstein and physical metrics should be 
more complicated than conformal; i.e., the so-called "disfor- 
mal transformation" in which certain directions are picked 
out for additional dilation or contraction (Bekenstein 1993, 
Bekenstein & Sanders 1994). An initial proposal for such a 
theory (Sanders 1997) invoked a non-dynamical vector field, 
with only a time component in the preferred cosmological 
frame, to provide this additional stretching. The disformal 
coupling was combined with an aquadratic Lagrangian for 
the scalar field to yield the MOND phenomenology. 

However, the non-dynamical aspect of the vector fleld 
violates general covariance making it impossible to deflne a 
conserved energy- momentum tensor (Lee, Lightman, & Ni 
1974). This problem led Bekenstein (2004) to construct a 
tensor-vector-scalar theory (TeVeS) with a fully dynamical 
vector fleld; this theory, while yielding MOND phenomenol- 
ogy in the weak field limit, is fully covariant, produces lens- 
ing at the same level as GR with dark matter, and possesses 
no obvious anomalies (propagation of scalar waves is causal) . 
In the same vein, I proposed a bi-scalar tensor- vector theory 



2 R.H. Sanders 



(BSTV) in order to provide a cosmological origin of Oo and 
cosmological dark matter in the form of scalar field oscilla- 
tions with wave-length sufficiently long to prevent clustering 
on the scale of galaxies (Sanders 2005). 

Thus it appears that any viable theory of MOND as 
modified gravity will require fields in addition to the tensor 
field of GR- a scalar field to yield the MOND phenomenol- 
ogy (as a fifth force) and a vector field to facilitate the non- 
conformal coupling and adequate gravitational lensing. In- 
deed, Soussa & Woodard (2004) have provided an elegant 
no-go argument to the effect that no single metric-based the- 
ory yielding MOND phenomenology in the weak field limit 
can produce the necessary degree of gravitational lensing. 
The only other possibility for MOND as modified gravity is, 
then, a multi-field theory such as TeVeS. 

The purpose of the present paper is to consider the 
constraints imposed upon multi-field theories of modified 
dynamics by solar system phenomenology. The precession 
of the orbits of Mercury and Icarus, as well as limits on 
the variation of Kepler's constant, GMq, between the earth 
and outer planets implies that the total force law within the 
orbit of Neptune is inverse square to high precision, apart 
from those post-Newtonian corrections introduced by GR. 
This suggests that any fifth force in the inner solar system, 
in addition to preserving the WEP, is also precisely inverse 
square. Moreover, the absence of detectable post-Newtonian 
effects due to a scalar field tied to a cosmic rest frame, i.e., 
ether-drift effects, probably constrains the magnitude of the 
fifth force to be less than 10~* that of the normal gravity 
force. But, in the context of MOND, the anomalous force in 
the Galaxy, at the neighbourhood of the sun, would have to 
be comparable to the gravity force. The transition from a 
weak inverse square attraction in the inner solar system to a 
significant anomalous attraction at several thousand astro- 
nomical units (au) would seem to require the appearance of 
a non-inverse square acceleration in the outer solar system. 

This is interesting in view of the fact that a deviation 
from inverse square attraction beyond 20 au is suggested on 
the basis of Doppler data from the two Pioneer spacecrafts, 
the Pioneer anomaly (Anderson et al. 1998, 2001). The mag- 
nitude of this apparently constant anomalous acceleration 
(~ 8 X 10^*^ crri/s'^) is tantalisingly close to, although sig- 
nificantly larger than, the MOND acceleration (Turyshev, 
Nieto & Anderson 2005). I show here that an anomalous 
acceleration is an expected, and indeed predicted, aspect of 
multi-field theories of modified dynamics. This non-inverse 
square acceleration appears in the outer solar system and 
need not be, but can be, as large as the observed Pioneer 
acceleration. Although the discussion is general, I illustrate 
this by considering current examples of multi-field theories 
of MOND. 

TeVeS, with Bekenstein's initial trial free function, pre- 
dicts a deviation which is too large to be consistent with 
both the reported constraints on A(GAdQ) and the proba- 
ble limits on preferred frame effects. These contradictions 
are not fatal because the free function of the theory can be 
modified to produce an anomalous force consistent with the 
planetary and preferred frame constraints. Indeed the form 
of the free function required is also consistent with that de- 
manded by observations of extended galaxy rotation curves 
which are flat beyond the visible disk. In this case the pre- 



dicted anomalous acceleration appears beyond 100 au and 
is roughly oo/3. 

The biscalar tensor vector theory (BSTV) is a modifi- 
cation of TeVeS constructed, in part, to be consistent with 
the constraints on deviations from inverse square attraction 
and with the non-detection of preferred frame effects near 
the earth. It also predicts a constant anomalous acceleration 
beyond Uranus that depends upon the value of the scalar 
coupling strength. For values of the scalar field coupling con- 
stant below a critical value then the constant acceleration 
is also w ao/3 as in TeVeS, but for larger couplings, the 
constant acceleration can be significantly larger than ao and 
extend within the orbit of Neptune; i.e., the theory may be 
tuned to be consistent with the Pioneer effect. If so, however, 
it is then inconsistent with the reported limits on deviations 
fi'om 1/r^ attraction out to the orbit of Neptune. This is 
unavoidable because if the Pioneer effect is really present 
within 30 au, then it would be inconsistent with limits on 
variation of GMq between the orbits of the inner planets 
and the orbit of Uranus and Neptune- limits derived from 
spacecraft ranging to these two outer planets. Either this 
constraint on deviations from l/r^ in the outer solar sys- 
tem is too stringent, which is possible (Section 4.3), or the 
reported Pioneer anomaly has a standard explanation (not 
involving fundamental physics). A more radical possibility 
is that the Pioneer effect, and hence MOND, is not due to a 
modification of gravity but of the particle action (Milgrom 
1994). 



2 MULTI-FIELD THEORIES OF MODIFIED 

DYNAMICS 

2.1 General properties of multi-field theories 

I have emphasised that, in scalar-tensor theories of MOND, 
the relation between the physical and gravitational metrics 
cannot be conformal. This condition requires the introduc- 
tion of a vector field. A" , that points in the time direction 
in the preferred cosmological frame. If the physical metric 
Qfiu is related to the gravitational metric gr^i. as 

Snii^ = e~'^'^'^g^iu - 2sinh{2ri(l))A^A^. (1) 

then it may be shown that the scalar field enhances the de- 
flection of photons about a visible astronomical system ex- 
actly as it would be by appropriately added dark matter in 
the context of pure GR; i.e., relativistic and non-relativistic 
particles feel the same total weaJj field force. (Sanders 1997, 
Bekenstein 2004). Here 77 is a parameter describing the 
strength of the scalar coupling to matter and is related to 
the parameter k in Bekenstein's notation {rj'^ — k/Aiv). 

It is useful to discuss scalar-tensor theories of modified 
dynamics in the context of the Einstein frame where the 
scalar field may be considered to mediate a force, fs, in 
addition to the usual gravity force connected to the grav- 
itational tensor, the Einstein- Newton force /n', that is to 
say, in this frame, particle motion is generally non-geodesic. 
In such theories the phenomenology associated with MOND 
results from this "fifth force" which is, in the extragalactic 
domain, a non-inverse square force that dominates in the 
regime of low field gradients {fs = rjc^Vtl) < ao). 

The aspect of non-inverse square attraction requires a 
departure from the standard Lagrangian for scalar-tensor 
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theories (Ls = 4>'°'4',o) either in the form of the aquadratic 
theory with a non-standard scalar Lagrangian {F{Ls)) or 
a biscalar theory where one field couples to matter and the 
second determines the strength of that coupling (as in PCG). 
Each of these proscriptions may bo designed to provide to a 
scalar force about a point mass of the form 



VGMao 



(2) 



at least in the regime where fs < ao- The total weak field 
force would then be given by ft = fs + fN where /jv = 
GM/r^ is the usual Newtonian force; clearly fs given by eq. 
2 will dominate at accelerations below ao- 

TeVcS is an aquadratic theory in disguise, with a scalar 
field action that may be written as 



Ls = ii'nV^-Vct> + V(n). 



(3) 



This is the weak coupling limit of PCG (77 << 1), the 
AQUAL limit, where one may show that the kinetic term for 
/X vanishes. As written here, /i is an auxiliary non-dynamical 
field and is algebraically related to (V(/))^ via the potential 
function V{n). This relation, expressed in terms of the scalar 
force fs = TjV^x? is given by 



(4) 



where V' = dV/d/j, and Im = c^/ao- Here I will refer to 

-lmV'{jj,) as the free function of the theory, although this 
function has, as its basis in the Lagrangian, the potential 
of a possibly dynamical field, V(fi). In the weak field static 
limit eq. 3 leads to the well known Bekenstein-Milgrom field 
equation 



V- [/x(|/,|/ao)/s] =47rGp. 



(5) 



The function fi{x) as it appears in cq. 5 docs not have 
the same meaning as fl in the original MOND prescrip- 
tion (/tAt(|/t|/ao) = /iv) or in the single-field Bekenstein- 
Milgrom non-relativistic theory (Bekenstein 2004). Eq. 5 
applies only to the scalar component of the force. In the 
context of such multi-field theories, not all forms of fi arc 
realisable from sensible single- valued forms of /i (Zhao & 
Famaey 2005). 

The most obvious, and simplest, choice for the free func- 
tion would be 



Voo/ 



(6) 



or ij,{x) = X for all x. This corresponds to Vi^ii) = (3Z^) 
and leads to a scalar force of the form of eq. 2 at all r; of 
course, the Newtonian force dominates for ft > ao. The 
rotation curves of spiral galaxies would be asymptotically 
flat as in MOND and would satisfy a mass-rotation velocity 
relation (TuUy-Fisher) of the form w"* oc M. We see below, 
however, that such a theory is inconsistent with the observed 
form of galajcy rotation curves as well as tight constraints on 
deviations from inverse square attraction in the inner solar 
system. 
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Figure 1. The log of the Newtonian force, /jv and scalar force, 
fs, in units of 10~* cm/s^ plotted against the log of the radial 
distance from a point mass in units of the MOND radius (rm = 

^GM/a„). The solid curve is the Newtonian force. The dotted 
curve is the scalar force for TeVeS with the free function originally 
taken by Bekenstein (2005). The siiort dashed curve is the scalar 
force resulting from the free function suggested by Zhao &; Famaey 
(2006) and the long dashed curve is the scalar force corresponding 
to eq. 8 here. In all cases r] = 0.01 



2.2 Rotation curve constraints on fifth force 

theories 

A more complicated scalar field Lagrangian is provided by 
the Bekenstein free function 



/s \ ^ 1 2 r 2 



(7) 



(here n as defined by eq. 5 differs by a factor of 7]^ from 
Bekenstein's definition). This yields a scalar force illustrated 
by the dotted curve in Fig. 1 where we see a return to 
attraction at high accelerations (here 77 = 0.01). For fs/ao < 
10~* this is equivalent to scalar force dependence provided 
by eq. 6. 

This free function, as well as that described by eq. 6, 

is unacceptable in that the form of the observed rotation 
curves of spiral galaxies implies that the scalar force cannot 
continue to increase smoothly as 1/r for accelerations near 
ao; the resulting rotation curves decline too slowly to the 
asymptotically constant value. This has been demonstrated 
for the Milky Way galaxy and for the well-studied spiral 
galaxy, NGC 3198 (Famaey & Binney 2005), and it is gen- 
erally true (Zhao & Famaey 2005). Fig. 2 shows the Newto- 
nian rotation curve (solid curve) resulting from a spherically 
symmetric meiss distribution of galaxy scale mass (lO^Afcs), 
an exponential sphere with a length scale of 2 kpc. The dot- 
ted curve is the rotation curve resulting from TeVeS with 
the free function described by eq. 7. The slow decline to the 
asymptotic value is evident. 
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Figure 2. Rotation curves for a spherical galaxy resulting from 
TeVeS with two alternative forms of the free function. The density 
distribution is taJsen to be exponential with a scale length of 2 
kpc and the total mass is 10^ ^Mq. The dashed line corresponds 
to the original Bekenstein free function (eq. 7) and the dotted 
curve to that of the free function considered here (eq. 8). The 
solid curve is the Newtonian curve. 



A free function given by 



(^) = [m' + 2/][l + V]-'[l - 2Zn(l - 



(8) 



would give rise to a radial dependence for the scalar force 
about a point mass of the form shown by the long dashed 
curve shown in Fig. 1. This is qualitatively similar to that 
originally suggested for the aquadratic stratified theory 
(Sanders 1997) and to the radial force dependence provided 
by the free function proposed by Zhao & Famaey (2006) 
(short dashed curve). Here we see that for total acceler- 
ations greater than oo the scalar force becomes constant, 
/s w Op « oo/3 before resuming 1/r^ dependence at larger 
total accelerations. The MOND interpolating function, fx 
(for spherical symmetry) corresponding to the free functions 
of Bekenstein, Zhao-Famaey, and eq. 8 are shown in Fig. 3. 
Consistency with observed galaxy rotation curves requires a 
transition to the Newtonian regime at least as rapid as that 
provided by the Zhao-Famaey free function. 

The relative merit of various free functions is not the 
topic here; the point is that radial dependence of scalar force 
must become rather flat at accelerations larger than oo to 
be consistent with galaxy rotation curves. This is evident in 
Fig. 2 where the rotation curve resulting from TeVeS with 
the free function described by eq. 8 is shown by the long 
dashed curve; this is more consistent with observed rotation 
curves which are generally flat beyond the visible disk. As 
will be shown below, this same qualitative behaviour is also 
required to meet the constraints on deviations from inverse 
square attraction in the inner solar system while avoiding 
observable preferred frame effects. 



Figure 3. The interpolating function fi, in the MOND prescrip- 
tion (/t/i(|/t|/ao) = /at) corresponding to the Bekenstein free 
function (eq. 7, dotted curve), the Zhao-Famaey free function 
(dashed curve), and the free function considered here (eq. 8, solid 
curve) . The indicated form is only strictly valid for spherical sym- 
metry, and in all cases r/ = 0.01 As emphasised by Zhao & Famaey 
(2006) with the Bekenstein free function the transition between 
the MOND and Newtonian regimes is too gradual to yield agree- 
ment with observed rotation curves. 



The BSTV theory has been designed to produce a radial 

dependence of the scalar force similar in form to that given 
by TeVeS with the modified free function, eq. 8, (dashed 
curve in Fig. 1), and therefore, rotation curves of the ob- 
served form (see Fig. 1 and 2 in Sanders 2005). There is, 
however, one important difference: Because the field deter- 
mining the strength of the scalar coupling in this theory, 
equivalent to ji, is dynamical (unlike TeVeS where it is an 
auxiliary field only) the value of the constant scalar accelera- 
tion near fs = ao depends not only upon the scalax coupling 
strength, r;, but also upon the value of the source mass and 
its distribution. This can make an important difference in 
outer solar system phenomenology. 



3 SOLAR SYSTEM CONSTRAINTS 
3.1 Planetary motion 

The most sensitive natural gravity probe in the inner so- 
lar system is provided by the orbit of Mercury. As is 
well-known General Relativity found its first experimental 
success in providing a non-Newtonian explanation of the 
anomalous precession of this orbit. The precise prediction is 
A9 = 43.03" per century, and the present observational re- 
sult agrees with this to better than 0.05" per century (Will 
2001). That is to say, precession resulting from any addi- 
tional non-Newtonian effect must be less than this limit. 
This provides a strong constraint on long range fifth force 
models. 



Solar system constraints on MOND 5 



Here I will parameteriso an additional non-Newtonian 
force in terms of a constant acceleration Op. It is straightfor- 
ward to demonstrate that the precession introduced by such 
a constant acceleration would be given by 



Table 1. Planetary constraints on a constant anomalous acceler- 
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GMp 



(9) 



where e is the eccentricity of the orbit and a is the semi- 
major axis. This may be rewritten as 

= 6.5(1 - e^)^ (^^zg^) (^) arcsec/century 

(10) 

where V is the mean orbital velocity. If, for Mercury, Qp < 
0.05" /century, eq. 10 would imply that Op < 4.0 x 10"^° 
cm/s^. That is to say, any constant anomalous acceleration 
present at the distance of Mercury from the sun, must be 
200 times smaller than the reported Pioneer acceleration. 

The precision of measured planetary precession de- 
grades rapidly for the other terrestrial planets, but the as- 
teroid, Icarus, remains a useful probe because of its near 
earth passage (a = 1.08 au) and its high eccentricity 
(e=0.83). Here, the precession predicted from GR is 10.3 arc 
sec/century and the observed precession is Qp = 9.8 ± 0.8 
arc sec/century (Weinberg 1972). Taking 0.8 arc sec/century 
as the limit on precession due to a constant acceleration we 
find, from eq. 10 that Up < 6.3 x 10~* cm/s^. 

Beyond Icarus, the tightest constraints on deviations 
from attraction are provided by limits on the varia- 
tion of Kepler's constant, Kp = GMq. If a variation, AKp 
is detected between two planetary orbits, and if this is pa- 
rameterised by the presence of a constant acceleration, Op, 
then 



' / 2 2-,-! 

-Kp(r2 — ri ) 



(11) 



where ri and r2 arc the distances from the sun of the closer 
and more distant planets respectively (assumed to be on 
circular orbits). An additional scalar force described by eq. 
2 (in addition to the Newtonian force) would result in a 
variation of Kepler's constant given by 



{rm = y/GMs/ao « 7700 au ). 

Since the advent of interplanetary spacecrafts, the dis- 
tances to the planets are known to high accuracy. This pro- 
vides strict limits on the variation of Kepler's constant be- 
tween the earth and the planet in question, i.e.. 



lAP 3 Ar 
2~P~ 2~ 



(13) 



where AP/P is the uncertainty in the period and Ar/r is 
the uncertainty in the solar distance. For example, from the 
Viking mission to Mars it is known that the uncertainty in 
the difference in the orbital radii of Earth and Mars is less 
than 100 m. Moreover, the difference in the orbital periods 
between the Earth and Mars is known to better than 7 parts 
in 10". By eq. 13 this implies that AKp/Kp < 2 x 10"^ 
and hence, from cq. 11, ap < 0.1 x 10^** for any constant 
acceleration present between the orbits of the Earth and 
Mars (Anderson et al. 2002). 

Similar observational limits resulting from Pioneer and 



atlon 


Object 


Distance (AU) 


Method 


Op (10-® cm/s^) 


Mercury 


0.39 




0.04 


Icarus 


1.08 




6.3 


Mars 


1.52 


^KplKp 


0.1 


Jupiter 


5.2 


/^Kp/Kp 


0.12 


Uranus 


19.2 


AKp/Kp 


0.08* 


Neptune 


30.1 


I^KpjKp 


0.13* 



The final column is the upper limit on a constant anomalous 
acceleration determined from planetary orbits via the indicated 
method (Hp, planetary precession; AKp/ Kp, variation of Kepler's 
constant). The constraints imposed by the orbits of Uranus and 
Neptune (marked with asterisks) are controversial (sec text), but 
if valid would be inconsistent with the Pioneer anomaly as a mod- 
ification of gravity. The distance given is the semi-major axis of 
the orbit.. 

Voyager flybys constrain AKp/Kp between the inner plan- 
ets and Jupiter, Uranus and Neptune to be less than .12, 
0.5, 2.0 times lO"'' respectively (Anderson et al. 1995). The 
corresponding limits on a constant acceleration are .26, .08, 
and .13 in units of 10"* era/ . We see that these limits for 
Urarms and Neptune are inconsistent with the reported Pio- 
neer anomaly if the anomaly is present at distances beyond 
20 au. These results are summarised in Table 1. 

There is disagreement over these outer solar system con- 
straints (Section 4.3), but in any case, it is clear that the 
1/r dependence of a fifth force cannot continue into the in- 
ner solar system- certainly not to within the orbit of Mars 
{T/vm. = 2 x 10"'')- because here the total gravitational field 
is so nearly inverse square. A force law of the form of eq. 2 
would resuh in AKp/Kp = 7 Ax IQ'^ (eq. 12) or 30000 times 
larger than the observed limit. This suggests that in TeVeS 
the radial dependence of scalar force must be quite precisely 
1/r^ certainly within the orbit of Mars. As we sec from Fig. 
1 this is consistent with the free function given by cq. 8 or 
by Zhao & Famaey- forms which arc also consistent with 
observed galaxy rotation curves. Note that a theory can be 
constructed in which the net scalar force vanishes within oo. 
This is highly contrived (involving two scalar components, 
one attractive and one repulsive), but implies that, in all 
that follows, one should add the condition that the scalar 
force is a monotonically decreasing function of radius. 



3.2 Post-Newtonian constraints 

Given that the total force must be quite precisely 1/r^ in 
the inner solar system, it is reasonable to suppose that the 
scalar force, in the high acceleration regime, is also 1/r^ as 
in Brans-Dicke theory. Then one may ask if there is any re- 
striction on the ratio of the weak field scalar to Newtonian 
forces in the inner solar solar system /s//jv. In Brans-Dicke 
theory, /s//jv — rf = l/(2a; -|- 3) (w is the Brans-Dicke 
measure of scalar coupling strength) . Because both fa and 
/jv are inverse square in the weak field limit, there is no 
restriction on the ratio of forces, or to, from weak field phe- 
nomenology; the restrictions appear at the post-Newtonian 
level. 

In isotropic co-ordinates the metric about a point mass 
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may be written as 



1+27 



GM 



dr 



(14) 

where the coefficients 7 and {3- the Eddington-Robertson 
parameters- describe the lowest order relativistic deviations 
from Newtonian inverse square gravity (post-Newtonian). In 
GR 7 = B = 1 precisely, and in Brans-Dicke theory it may 
be shown that /3 = 1. In fact, this is true of any conformally 
coupled scalar-tensor theory, = V'(0)Smi' provided that 



<(0) = [^^'(0)]= 



(15) 



(see Appendix). But, as noted in the Introduction, because 
of the conformal relation between the physical and gravita- 
tional metrics, there is no enhanced deflection of photons 
due to the scalar field, while non-relativistic particles do re- 
spond to an enhanced force. This is reflected in the fact 
that the post-Newtonian parameter 7 = (a'-|-l)/(a;-|-2) ^ 1 
in Brams-Dicke theory. In general, 7 7^^: 1 in in conformally 
coupled scalar-tensor theories. 

A disformal transformation of the form of eq. 1 has 
been discussed by Giannios (2005) in the context of TeVeS. 
It is equivalent to multiplying different components of gij,i, 
by separate functions of (j>; i.e., gtt = ^l){<j})gtt and Qrr = 
x{(j>)grr- In such theories, it is the case that 7 = 1 if 



x'(o) = -v'(o) 



(16) 



(see Appendix). If both conditions 15 and 16 are met and 
if ^4*^ = (i.e., the vector does not develop a non-zero 
radial component), then /3 = 7 = 1. It is easy to ver- 
ify that the particular transformation provided by eq. 1, 
where ipi^cj)) = exp{2'q(f>) and x(<^) — e2;p(— 2r7(/)), satis- 
fles conditions 15 and 16. Therefore, for any tensor-vector- 
scalar theory in which the gravitational and physical met- 
rics are related according to eq. 1, there is no restriction on 
fs/ Jn at post-Newtonian level as described by the standard 
Eddington-Robertson parameters. This is true of the classi- 
cal stratifled theories (Ni 1972), of the stratified aquadratic 
theory (Sanders 1997) and of TeVeS assuming A"" = 0. These 
theories are consistent with a wide range of observed phe- 
nomena from deflection of starlight by the sun to radar echo 
delay. It is the gravitational preferred frame effects that are 
threatening for such theories. 



3.3 Preferred frame constraints 

Multi-field theories of MOND must contain a normalised 

cosmic vector field to provide the disformal transformation. 
The direction of the vector is determined primarily by the 
universal mass distribution and, in a FRW metric, points 
in the positive time direction. The equations of motion in 
a gravitational field take their simplest form in the cosmic 
frame where only the time component of the vector field is 
non-zero. For a frame in relative motion, such as the solar 
system, space components of the vector field develop non- 
zero values, and this aifects the motion of particles; i.e., ether 
drift effects must appear at some level; i.e., such theories 
violate the Lorentz invarianee of gravitational dynamics. 

Post-Newtonian preferred frame effects in conservative 
theories are quantified by two parameters, ai and 02 (Will 
& Nordtvedt 1972). These modify the effective Lagrangian 



that describes the gravitational dynamics of N-body systems 
by adding terms such as 



and 
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(18) 



where Va is the velocity of particle A with respect to the 
preferred cosmic frame, w is the velocity of the inertial frame 
with respect to the cosmic frame, and tab is the unit vec- 
tor along rAB- A non-zero value of ai would lead to ef- 
fects such as a polarisation of the earth-moon orbit and 
is constrained to be less than 10"'* by Lunar Laser Rang- 
ing (Miiller, Nordtvedt & Vokrouhlicky 1996). The 0:2 term 
quantifies effects such as periodic variation in the effective 
gravitational constant (with twice the orbital or rotational 
frequency of the system) or an ether drift torque acting on 
a spinning body. This is constrained to be less than 10~^ by 
the near-alignment of the sun's rotational axis with that of 
the solar system (Nordtvedt 1987). 

Calculation of the predicted values of ai and 02 must 
be done for each particular theory. Here to keep the dis- 
cussion as general as possible, I provide estimates of the 
preferred frame parameters in tensor-vector-scalar theories 
by heuristic arguments. 

In the historical Lagrangian-based stratified theories 
such as that of Ni (1972) there is one dynamical field, (f>, 
a noii-dynamieal tensor (Minkowski) describing the back- 
ground geometry, and a non-dynamical vector field; that is 
to say, the gravitational force is supposed to be mediated 
only by the scalar field disformally coupled to an a priori 
geometry described by the Minkowski metric. Here it may 
be shown that a2 = but ai = —8 in sharp contradic- 
tion with the LLR result, not to mention earlier constraints 
on the diurnal and annual variation of the the gravitational 
constant. 

In the predecessor to TeVeS, the aquadratic stratified 

theory (Sanders 1997), there are two dynamical fields, a 
scalar and the Einstein metric, in addition to a non-dynamic 
vector field. That is to say, the gravitational force is medi- 
ated not only by the scalar field, but also by the Einstein 
metric which is locally insensitive to motion through the 
cosmic frame. Here the preferred frame effects would ap- 
pear through the contribution of the scalar to the physical 
metric (ala eq. 1). In the limit where the scalar coupling, 
rj, vanishes, the theory reduces to GR and in GR there are 
no preferred frame effects. Therefore observational limits on 
preferred frame effects must place an upper limit upon rj. 

This can be made more definite by noting that the equa- 
tion of motion for the scalar field, in the high acceleration 
limit where the Lagrangian is standard, has the form 



<;/3i 



■.13 



: 47rG77r^.[5'"' + {1 + 6-^""*)^^^" 



(19) 



where ff^, is the usual energy momentum tensor in the phys- 
ical frame (Bekenstein 2005). In a frame moving with veloc- 
ity w with respect to the cosmic frame, the source, to order 
ui^/c^ would take the form 47rGr;p[l -|- {w^ +w ■ t;)/2c^], but 
Qfiv would contain no such terms to post-Newtonian order. 
Since fs = riVcj), I conjecture that ai is suppressed, relative 
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to its value in pure stratified theories, by « fs/fN- If so, 
this would constrain fs/fN < 10~* 

Similar arguments would apply in a theory such as 
TeVeS where all three fields, including the vector, are dy- 
namical; that is, wc would expect post-Newtonian preferred 
frame effects to project into the solar system via the scalar 
field which is tied to the cosmological frame. There is, how- 
ever, an additional effect because the vector field contributes 
directly to the source of the Einstein tensor, G^^. This con- 
tribution is not negligible because of the presence, in the the- 
ory, of a Lagrangian multiplier function included to enforce 
a normalisation condition on the vector field, Aij,A^ = — 1. 
The additional term in the energy-momentum tensor is then 
— XA^Au which, from the vector field equation becomes 
« AtiKGn f>Ai_iA,j where K is & new parameter associated 
with the vector field (a coupling strength parameter), and 
Giv is the locally measured gravitational constant (the pri- 
mary effect of the vector field in the weak field limit is a 
rescaling of the gravitational constant with respect to its 
cosmological value, G, i.e., Gn ~ G{1 — K/2)~^). Therefore, 
given a Lorentz transformation of the cosmic vector field to 
a moving frame we see that the source of G^i. contains terms 
proportional to KGn pw^ /<? ■ This would constrain K to bo 
less than qi or a2 (say K < 10"^) but would have no pro- 
found effect on weak field phenomenology. 

In summary then, we can say that the suppression of 
the likely preferred frame effects such as polarisation of the 
earth-moon orbit will probably require that fs/fN < 10^* 
This is comparable to the current reported constraints on a 
scalar force in the context of Brans-Dicke theory (Bertotti 
et al. 2003). But I re-emphasise, this is a heuristic argument 
and a proper calculation should be done. 



4 CONFRONTATION OF MULTI-FIELD 
THEORIES WITH SOLAR SYSTEM 
PHENOMENOLOGY 

Summarising the above discussion, we have seen, first of 
all, that the total weak field gravitational force in the solar 
system is inverse square to high precision, at least within 
the orbit of Mars. This implies that any component of the 
gravitational force, in addition to the Einstein-Newton force, 
should also be precisely inverse square. At the same time, 
post-Newtonian preferred frame effects would seem to re- 
quire that any additional inverse square force due to a scalar 
field should be smaller than 10~* of the Einstein-Newton 
force. It is not trivial for a theory to satisfy these two con- 
straints. 

4.1 TeVeS 

In TeVeS, the relevant parameter, which determines the 
strength of the scalar force is r;; i.e., fs/fN ~ rf in the 
limit where fs » ao. The preferred frame considerations 
would then require the r; < 10~^. In Fig. 4 the dotted curve 
shows the anomalous force (the non component of the 
total force) resulting from Bekenstein's initial free function 
(cq. 7), with r/ = 0.01, compared to the total force within 
the solar system (dashed curve). The points are the limits 
on a constant anomalous acceleration from planetary mo- 
tion discussed in section 2. The solid bar represents the Pi- 
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Figure 4. The dashed curve is the log of the the total force 
(/( = fs + /jv)) in units of 10~® cm/s-^ plotted against the log of 
the radial distance from the sun in astronomical units for TeVeS. 
The dotted curve is the anomalous force (the non-inverse square 
force) for Bekenstein's initial choice of free function with rj = 0.01 
(cq. 7). The long dashed curve is the same but with t] = 0.1. The 
solid curve is the anomalous force resulting from the modified 
free function (eq. 8). The points show the limits on a peculiar 
acceleration implied by planetary motion as discussed in Section 
2, and the solid bar is the Pioneer acceleration. 

oneer anomalous acceleration. Here it is obvious that the 
theory, with this choice of free function and scalar coupling 
strength, strongly violates these limits on deviations from 
inverse square attraction well into the inner solar system. 

The deviation from inverse square attraction is less se- 
vere if the scalar coupling constant, rj is larger. The anoma- 
lous force resulting when 77 = 0.1 is shown by the long dashed 
curve; this would appear to be roughly consistent with the 
constraints on deviations from inverse square within the in- 
ner solar system. But then the scalar force is only 0.01 of the 
Newton-Einstein force, and the theory would probably evi- 
dence local preferred frame effects at least 100 times larger 
than the present limits. 

The modified form of the free function (eq. 8) with 77 — 
0.01, gives an anomalous force shown by the solid curve in 
Fig. 4. It is obvious that this is consistent both with the 
planetary constraints on deviations from 1 /r^ attraction out 
to Neptune and with the avoidance of local preferred frame 
effects required by fs/fN < 10^'*, but a constant anomalous 
acceleration « ao /3 does appear beyond 100 au. 

4.2 BSTV 

BSTV is in part designed to satisfy these solar system con- 
straints but at the expense of adding a new parameter e> r) 
(the parameter is not necessary an additional; it may be 
identified with the vector coupling strength) . Here there are 
two explicitly dynamical scalar fields- one, 0, that couples 
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to matter and the second q which determines the strength 
of the coupling. In terms of the scalar field gradient, the 
quasi-static field equation is 



(20) 



where in the high acceleration limit; i.e., q saturates 

at a small value in this limit (note that /U = /2rf with as 
defined in eqs. 3 and 5). Given that r) is the strength of the 
scalar field coupling (as in cq. 1) and that the scalar force is 
fs = ric^\7(f>, the theory is designed to yield a precisely 
force in the inner solar system with fs/ Jn 2rf /e^ in the 
limit where fs » ao- Thus, for this theory, the avoidance 
of preferred frame effects near the earth would require that 

But there is another significant difference with TeVeS. 
The relation between the coupling strength field, q, and the 
scalar force is no longer algebraic (as in eqs. 7 and 8) but is 
differential and given by 



V q 



■qV<P-V^ = V:{q) 



(21) 



where Vs (q) is now an effective potential involving the cos- 
mic time derivative of the scalar field. 

It is instructive to view this equation in uni t-less for m 
by defining y = q/r] and x = r/vm with r-m = •s/GMJoo = 
^/tsIm (ts ~ 2GM/(? is the Schwarzschild radius). Given 
that V {q) ^ 2g^c/>^/e^ in this regime (q « r;), that ao = 
^/V2ri\(j>\/t, and that V(t> = 'qGM /{(?r'^q'^), then, in the case 
of spherical symmetry, eq. 21 becomes 



2 'm _2 
T) — V y = 

Vs 



x^y^ 12^' 



(22) 



From this it is obvious that, in the limit of weak coupling, 
where iflu/rs « 1, the relationship between q'^ and V0 
is effectively algebraic as in TeVeS (the theory approaches 
its AQUAL limit). In this case the radial dependence of 
the scalar force is similar to that shown in Fig. 1 (dashed 
curves). This is consistent with galaxy rotation curves and 
would produce an anomalous acceleration in the solar sys- 
tem similar to that of Fig. 4 (solid curve); i.e., it satisfies all 
planetary constraints on deviations from 1/r^ attraction. 

For a given value of r] the condition for weak coupling 
provides a lower limit on the source mass. Whenever 



M > M^^ 6r;- 



G 



5 X 10^^r;^Me 



(23) 



then the weak coupling limit applies (here Ih = /Ho is 
the Hubble radius and I have taken Im = 61h). If r? w 2 x 

lO"^'^ this critical mass would correspond to a few solar 
masses. In other words, for larger mass, the weak coupling 
limit applies, and the form of the scalar force (as a function 
of r/rm) is frozen as in TeVeS. But for smaller masses, the 
full differential equation (eq. 22) must be solved and the 
solution depends upon the source mass (the presence of a 
critical mass in PCG was pointed out by Bckcnstcin 1988). 

This is illustrated in Fig. 5 where we sec the scalar force 
as a function of the scaled radius {r/rm) for r? = 2 x 10"^^ 
and for various values of the source mass. The object is 
placed in the galaxy acceleration field near the position of 
the sun where g « (the scalar force is comparable to the 
Newtonian force at large distance from the star). Here, if 
M > 5OM0, the solution is fixed at the weak-coupling limit. 
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Figure 5. The log of the scalar force {fs, solid curves) and the 
total force {ft = fs + /iVi dotted curve), in units of 10~* cm/s^ 
plotted against the log of the radial distance from a point mass in 

units of the MOND radius rm, = GM/ao for the biscalar the- 
ory. The various curves correspond to the indicated values of the 
source mass. The curves converge for M > 50Mq corresponding 
to the weak coupling limit (the AQUAL limit) of the theory. 



For smaller values of the mass, the solution, and in particu- 
lar the value of the plateau acceleration, depends upon the 
source mass. If the source mass is IMq, then the plateau ac- 
celeration is 8 x 10^** cm/s^. That is to say, unlike the weak 
coupling or AQUAL limit, the constant anomalous acceler- 
ation in the outer solar system can be significantly larger 
than the MOND critical acceleration. 

With r; = 2 X 10-12 and = 2 x IQ^r]^ (Sanders 2005), 
the resulting anomalous force (solid curve) is compared to 
the total force (dashed curve) in Fig. 6. Here it is evident 
that, within 20 au the scalar force is precisely 1/r^ and 10-'' 
less than the Newtonian force. Moreover, the theory in this 
form produces a constant anomalous force that is consis- 
tent with the Pioneer acceleration but inconsistent with the 
the reported limits the variation of Kepler's constant out to 
Uranus and Neptune. Taking the parameters of the theory 
to be ?7 = 0.9 X 10-'^ and = 4 x lO''?;^ pushes the the- 
ory back to the AQUAL limit and produces the anomalous 
force shown as the dotted curve. This is consistent with all 
reported planetary constraints on inverse square attraction 
and preferred frame effects as is TeVeS with the revised free 
function. In fact, it is identical in this respect to TeVeS with 
the revised free function (eq. 8). 



4.3 The Pioneer anomaly 

TeVeS with the free function modified to be consistent 
with galaxy rotation curves (as in oq. 8) predicts a con- 
stant anomalous acceleration beyond 100 au with magni- 
tude w 3 X 10-® cm/s^. The same is true of BSTV in the 
limit of weak scalar coupling, r] < lO-'^^. Therefore the the- 
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Figure 6. The dashed curve is the log of the total force {ft = 
is + /n), in units of 10~* cm/s^ plotted against the log of the 
radial distance from the sun in astronomical units for the biscalar 
theory. The solid curve is the anomalous force (the non-inverse 



square component) for a scalar coupling r/ = 2 X 10 



-12 



(with 



e = lAlrj) and the dotted curve is for rj = 0.9 X 10^-'-^ (witli 
e = 200?}) . As in Fig. 4, the points show the planetary constraints 
on an anomalous non inverse square attraction. 



orics may be constructed to be consistent both with galaxy 
rotation curves and with planetary constraints on at- 
traction within the orbit of Neptune. In any case, an anoma- 
lous acceleration, dp, in the outer solar system is inevitable, 
provided that the scalar force is a monotonically decreasing 
function of radius. Consistency with galaxy rotation curves 
appears to require that ap w 0.3ao as a lower limit. For 
BSTV, however, because the field determining the efi^ective 
strength of the scalar coupling, q (or ^) is dynamical, it is 
possible that Op > ao beyond 20 au while, in the outskirts 
of galaxies ap w 0.3ao. The same would probably be true 
of TeVeS with a dynamical /i. It is therefore tempting to 
identify the predicted constant acceleration with the Pio- 
neer anomaly. 

This is problematic because, as we see in Fig. 5, the Pi- 
oneer anomaly itself is inconsistent with reported limits on 
the variation of Kepler's constant out to Uranus and Nep- 
tune (Anderson et al. 1995). However, these stated limits 
may be overly stringent because they are based only upon 
single spacecraft ranging measurements to these planets, and 
the formal uncertainties in the distances are almost certainly 
too optimistic. It should also bo kept in mind that both 
Uranus and Neptune have not completed a single orbit pe- 
riod since the advent of precise astronomical positioning in- 
strumentation, and, therefore, their orbits are poorly known 
(Standish 2004 and private communication 2005). 

More recently, it has been claimed that such a large 
anomalous acceleration, if present beyond 20 au, would lead 
to secular and short period signals in the orbits of the outer 
three planets- signals large enough to have been detected 



given the present levels of accuracy (lorio 2006) . The oppo- 
site conclusion has been reached by Page, Dixon and Walen 
(2005) who propose using distant asteroid orbits as a test 
for the Pioneer effect. It would seem fair to conclude that 
there is an lack of agreement about the nature of the grav- 
itational field (as probed by planetary orbits) in the outer 
solar system. 

This is an important issue. If planetary motion beyond 
20 au is inconsistent with the presence of the constant Pi- 
oneer acceleration, then the Pioneer anomaly is not due to 
a modification of gravity in the usual sense. If the plane- 
tary motion is consistent with the Pioneer anomaly, then it 
remains possible that this reported constant acceleration is 
due to the effect of a fifth force which becomes evident at 
low accelerations, as in relativistic theories of MOND. 



5 CONCLUSIONS 

Solar system phenomenology, in particular the tight limits 
on deviations from inverse square attraction and the ab- 
sence to high precision of local preferred frame effects, places 
strong constraints on multi-field theories of modified dynam- 
ics as modified gravity. Specifically, any fifth force mediated 
by a scalar field must also be inverse square to high preci- 
sion in the inner solar system, at least within the orbit of 
Mars, but smaller than about 10^* of the Newton-Einstein 
force to avoid producing observable preferred frame effects. 
I re-emphasise that this constraint upon the magnitude of 
the scalar force is only an estimate based upon heuristic 
arguments; a proper calculation of the preferred frame post- 
Newtonian parameters for TeVeS should be done. It is clear, 
however, that in these theories preferred frame effects, such 
as a polarisation of the lunar orbit, should appear at some 
level. 

In the Galaxy, at the position of the sun, the galactic 
gravitational acceleration is on the order of MOND acceler- 
ation oo ~ 10~® cm/s^ . In the context of multi-field theories 
of MOND this would imply that a fifth force acceleration 
is probably about as large as the Newtonian acceleration, 
or, in terms of scalar-tensor theory, /s//jv ~ 1. In other 
words, fs/fN must grow from 0.0001 within the orbit of 
Neptune to about one at a distance of 800 au where the 
galactic gravitational acceleration becomes comparable to 
the solar attraction. Therefore, going outward from the sun 
to the galactic environment, the scalar force must appear 
as an anomalous non-inverse square acceleration (provided 
that the scalar force dependence on radius is monotonic). 
This effect, first noted for the stratified aquadratic theory 
(Sanders 1997), is an inevitable consequence of multi-field 
theories and is evidenced both by TeVeS and BSTV. 

In order to meet the solar system constraints of pre- 
cise inverse square attraction and the absence of preferred 
frame effects the scalar force must have the form demon- 
strated by the dashed curves in Fig. 1; i.e., there must be a 
transition region between 1/r and attraction where the 
acceleration due to the scalar field is more-or-less constant. 
Therefore, not only must an anomalous acceleration appear 
in the outer solar system (certainly beyond 100 au) but it 
must also be, to lowest order, constant with radius between 
100 and 1000 au. 

Both BSTV and TeVeS with the modified free function 
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can provide this transition, but there is an important differ- 
ence. In BSTV, it may be shown that for vanishing coupUng 
strength {ri —> 0) the Laplacian of q, the coupling strength 
field, may be neglected in the field equation for q, and the 
theory becomes, in effect, an aquadratic theory as in ToVeS. 
However, for finite r/ the form of q and hence the scalar force, 
in general, depends upon the mass of the source and the cou- 
pling strength. For rj < 10~^° the AQUAL limit applies to 
galaxy scale masses but the full differential equation must 
bo solved for smaller masses. The practical consequence of 
this is that the plateau acceleration (where fs < ao) de- 
pends upon the source mass. For a galaxy scale mass this 
near constant acceleration can be ~ Oo/S (as required for 
rotation curves) but for a solar mass it may be near lOao 
if 77 « 10~^^. This possibility exists for TeVeS as well if the 
auxiliary field, n, is given its own dynamics by writing a ki- 
netic term proportional to lJi,aH'°' into the Lagrangian (this 
would provide a more familiar theory). Therefore it is pos- 
sible that the the predicted constant anomalous force could 
be identified with the Pioneer anomaly. 

But it is also evident from Figs. 4 and 6 that no the- 
ory of MOND as modified gravity can satisfy the reported 
limits on deviation from inverse square attraction at the 
orbits of Uranus and Neptune and be consistent with the 
Pioneer effect if the anomalous acceleration really does ap- 
pear at radii as small as 20 au. This is because the Pio- 
neer effect itself is inconsistent with these constraints. The 
constraints themselves arc controversial. Basically, the so- 
lar gravitational field in the outer solar system is not well- 
understood, and this calls for a reconsideration of the orbits 
of the outer planets in the presence of a non-inverse square 
acceleration. 

It should also be noted that Milgrom (1994) has pro- 
posed a basis for MOND as modified inertia in which the 
particle action is a non-local functional of the entire particle 
trajectory. This is a completely different approach from the 
modified gravity theories discussed here, and could account 
for the possibility that the Pioneer spacecrafts on hyper- 
bolic orbits feel the anomalous acceleration but the planets 
on more circular orbits do not. For this reason, it would be 
of considerable interest to determine if the Pioneer anomaly 
first appears at the point where a spacecraft is boosted from 
a bound to an unbound orbit. 

The significance of the Pioneer effect should not be un- 
derstated. It may constitute the first evidence on a scale 
smaller than galactic and extra-galactic that there is more 
to gravity than we have supposed. The question of whether 
or not the Pioneer acceleration is a new physical effect and, if 
so, where the anomalous acceleration first appears requires 
rearialysis of the existing Pioneer data and, on the longer 
term, new space missions to confirm (or not) this important 
result (Turyshev, Nieto & Anderson 2004). 

The possibility of such local tests is, in a sense, the 
"holy grail" of modified gravity theories (as is the direct de- 
tection of new particles for dark matter theories). In this 
regard, Bekenstein and Magucijo (2006) have demonstrated 
that, in the context of TeVeS, MOND tidal stresses become 
anomalous large near saddle points in the local solar system 
where the total weak field force approaches zero, between 
the earth and the sun, for example. Space missions with sen- 
sitive accelerometers might detect such effects. This would 
indeed be a spectacular confirmation, but non-detection at 



the predicted level would not be a falsification of general 
modified gravity theories for MOND. In the biscalar theory, 
for example, the coupling strength field, also being dynam- 
ical, does not respond immediately to changes in the scalar 
force; once in the deep Newtonian regime, one cannot re- 
turn to the MOND regime over relatively small distances. 
On the other hand, it does appear that any theory of MOND 
as modified gravity would require the presence of an anomar 
lous acceleration in the outer solar system (beyond 100 au) 
with a magnitude of at least a few tenths ao. 

I am grateful to Jacob Bekenstein, Moti Milgrom, 
and Dimitrios Giannios for very useful comments on this 
manuscript. 1 also thank Slava Turyshev for helpful remarks 
on the Pioneer anomaly and Miles Standish for comments 
on planetary constraints on the solar gravitational field. 



APPENDIX A: EDDINGTON-ROBERTSON 
POST-NEWTONIAN PARAMETERS IN 
THEORIES WITH DISFORMALLY RELATED 
METRICS 

Let us suppose that, in the preferred frame, the relation be- 
tween Qtt and gtt (the time-time components of the physical 
and Einstein metrics) can be written as 

9tt = '4'{4>)gu (Al) 

where tp is a general function of (p, the scalar field. Further 
take the Taylor expansion 

1, 



where a = ip'{0) and b = tp"{0). If the scalar field dynamics 
is described by the standard field Lagrangian, as it is for 
these theories in the inner solar system (i.e., Lg = 4>,a4>'" 
then it is the case that 

()> = -a (A3) 

r 

as is usual in scalar-tensor theories. Taking 

gtt = -1 + hit (AA) 
where, to second order 



ht. 



(A5) 



Then we find, to second order 



gtt = -l + 2{l + -)—-2{l + a' + —){ — ) . (A6) 
Redefining the mass M' = (1 -h %-)M, we then have 



gtt = l-2^-2{^Y[l±^^- 

By identification with eq. 14 we find 



1 + g" + V 
(1 + ^)^ ■ 



(A7) 



(AS) 



The condition for /3 = 1 is then h = a or 

V"(0) = V'(0)' (A9) 

This is obviously true if V'(0) = which would be the case 
for the particular transformation described by eq. 1. 
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For a conformal transformation it is also the case that 
ffrr ~ ip{4>)grr- Repeating the procedure above for grr to 
first order, we find 



1 - fflV2 



(AfO) 



(given that a? — 2/ (2a; + 3) this returns the usual Brans- 
Dicke result). 

Now suppose the transformation is disformal, and of a 
form which generalises eq. 1: 



{All) 



where «(</>) and w{(l>) are unspecified functions of the scalar 
field. To post-Newtonian order this is equivalent to multi- 
plying the time-time and space-space components of by 
different functions of the scalar field, e.g., gtt ~ ^{4')gtt and 
grr = x(4')grr where 



and 



^(<^) = u(<^) -f ™(0) 

x{4>) = u{<p) 



Taking 

X{^) = 1 + aV 
and repeating the above procedure for g, r we find 
_ 1 - aa'/2 



l + a'^/2 



Here 7 = 1 if a = a' or 



x'(o) = -^'(0). 



(A12) 
{A13) 
(AU) 

{A15) 

(AW) 



Rewriting conditions A9 and A16 in terms of the func- 
tions u and w (via eqs. f 2 and f3) we find that if 

[u'{0) + w'{0)f = u"{0) + w"{0) {A17) 



and 



2u'(0) = -w'{0) (A18) 



then /3 = 7 = 1 It is straightforward to confirm that the 
particular transformation given by eq. f satisfies these two 
conditions. 

It is important to note that Giannios (2005) has found 
two spherically symmetric static solutions for TeVeS: one in 
which A^ = everywhere (as assumed here), and a second in 
which the vector field develops a non zero-radial component. 
In this second case /3 7^ and depends upon the vector 
coupling strength parameter and the cosmological value of 
the scalar field. It is not yet determined if both solutions are 
stable. 
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